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We study jets of massless particles in A/" = 4 super Yang-Mills using the AdS/CFT correspondence 
both at zero and finite temperature. We set up an initial state corresponding to a highly energetic 
quark/anti-quark pair and follow its time evolution into two jets. At finite temperature the jets 
stop after traveling a finite distance, whereas at zero temperature they travel and spread forever. 
We map out the corresponding baryon number charge density and identify the generic late time 
behavior of the jets as well as features that depend crucially on the initial conditions. 
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I. INTRODUCTION 

In a hadronic collider such as the LHC, most interest- 
ing processes involve final states with jets of hadrons or 
equivalently, sprays of particles heading in roughly the 
same direction within a certain opening angle. The data 
available to define and reconstruct the jets are energy de- 
positions in the calorimeter and charged particle tracks. 
Ideally one can use this experimental data to identify 
an underlying hard scattering event that involved QCD 
partons in the weakly coupled high energy regimes or, 
alternatively, heavy standard model particles such as top 
quarks or newly produced non-standard model particles. 
Distinguishing the former from the latter is of course cru- 
cial to our ability to find and analyze new physics at the 
LHC. 

Despite its obvious importance, several aspects of jet 
physics are still in need of better understanding. A large 
portion of the current understanding of jets comes from 
perturbative techniques. For a recent review see Ref \X\. 
First, a jet has to be defined via one of several algorithms, 
all of which have their own benefits and problems. The 
properties of the jet then have to be compared to a per- 
turbative calculation of the hard scattering. The initial 
hard scattering however has to be first processed by a 
parton shower routine that evolves every hard parton into 
a jet of partons which eventually hadronize into baryons 
and mesons. The distribution of the latter can then be 
matched with the experimentally observed data. This 
compartmentalization of the theoretical part into three 
stages assumes that these processes, which are governed 
by physics at very different energy scales, happen inde- 
pendently. This is the statement of factorization. While 
resummed perturbation theory or soft coUinear effective 
theory give a theoretical handle on the showering, Monte 
Carlo simulations based on models (and not QCD) are 
often employed for this step. The hadronization is always 
handled by phenomenological models. One question we 
will address in this work is how the process of shower- 
ing occurs in a particular solvable toy model of strong 
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coupling dynamics. In particular, we study jets in large 
A^c, SU{Nc) TV = 4 super- Yang-Mills (SYM) at strong 
't Hooft coupling A where the process has a dual descrip- 
tion in terms of a string falling in a bd geometry. This 
may help us understand what aspects of showering are 
genuinely QCD and may also lead to better models of 
showering. 

A quite different experimental setting that prominently 
features jets are heavy ion colliders, such as the currently 
operating RHIC experiment as well as the heavy- ion runs 
that will take place part of the time at the LHC. The goal 
of these experiments is not so much to find physics be- 
yond the standard model but rather to explore the phase 
structure of QCD at high temperatures. At the temper- 
atures achieved in the RHIC collisions, quarks and glu- 
ons are thought to be liberated although not free. Their 
collective dynamics are well described by an almost per- 
fect fluid [5] [3]. Jets, which are produced by hard col- 
lisions within the expanding fireball, serve as a probe of 
the interior of the plasma. In this context one wants 
to know the rate at which jets lose energy and broaden 
when traversing the medium. These are strongly coupled 
non-equilibrium phenomena and are therefore challeng- 
ing for both the lattice and for diagrammatic approaches. 
Again, it is helpful to have a solvable toy model. We will 
study jets of massless quarks traversing a finite tempera- 
ture plasma in the same solvable toy model, M — A SYM 
at strong couphng. 

Our approach, both at zero and finite temperature, is 
to prepare at a state at time i = which is well described 
by a semi-classical wave- function. In the gravity dual, 
this corresponds to a classical string moving in the dual 
geometry. At i = we specify the initial configuration of 
the string and let it time evolve according to the equa- 
tions of motion that follow from the Nambu-Goto action. 
The time evolution of the string is dual to the evolution 
of baryon and energy densities in the corresponding field 
theory state. A similar approach has been used to study 
the properties of jets of massive matter in the same the- 
ory at zero temperature 0] [5] and finite temperature [^ . 
In this work, we compute the evolution of baryon density 
at both zero and finite temperature for the case of mass- 
less quarks. For concreteness, we consider states with 
two back to back jets. In particular, we consider string 
states where the string profile lies in one Minkowski spa- 
tial direction. However our analysis and results for the 



baryon density are easily generalized to the case where 
the string profile extends in more than one Minkowski 
spatial direction. 

For the classical description to be valid the string must 
carry large conserved charges. In particular, the energy 
of the string scales as vA. At finite temperature, this 
setup is very well suited for the question of studying the 
energy and momentum loss rates of an energetic projec- 
tile traversing the medium. Our calculation generalizes 
the calculation of the heavy quark energy loss rate of 
[SJ [7j to the case of massless quarks. The additional 
complication is that for massless quarks the plasma can 
transport the charge and so the initially localized charge 
distribution eventually diffuses; the string falls. While we 
don't have a conclusive answer on the loss rate and stop- 
ping distance yet, our simulations indicate that a recent 
analysis of this process based purely on the analysis of 
geodesies [8] is incomplete. One of the central results of 
our numerical investigations is that the endpoint motion 
indeed very quickly settles onto a lightlike geodesic for 
generic initial conditions. However the relation between 
the energy of the quark and the parameters specifying 
the geodesic can be more complicated than presented in 
0. 

At zero temperature our approach is new, even though 
the situation is quite similar in spirit to the case of 
finite temperature. At time i = we prepare a 
quark/ antiquark pair with enough kinetic energy to climb 
out of each others' potential well. We follow in detail 
the evolution of the baryon number density and can in 
principle do the same for energy densities. As each one 
of the energetic projectiles moves through the vacuum, 
both densities spread so that instead of two single well de- 
fined particles fiying back-to-back at the speed of light, 
we see narrow distributions of baryon number density 
within two cones around the directions of the original 
quark/ antiquark pair — a pair of jets. The dynamical 
mechanism underlying this jet formation from the field 
theory perspective seems to be the emission of a large 
number of soft gluons and quarks. This setup with an 
initial semi-classical state whose subsequent evolution is 
studied mirrors closely the usual factorization framework. 
An initial hard scattering event will have to produce the 
i = string configuration. We do not specify the dy- 
namics that gave rise to the initial configuration. In 
any case, this would have to be very different from QCD 
where the initial hard event is perturbative. The subse- 
quent showering is what we study using the gravity dual. 
The conformal theory itself doesn't hadronize, but grav- 
ity duals for confining gauge theories are known. In these 
cases only the deep interior of the geometry is modified. 
This will not modify the initial falling stage of the string 
motion. Factorization of the showering stage from the 
hadronization stage is guaranteed by bulk locality in our 
setup. 

This approach is rather different from the recent work 
in [3]. There the authors study the evolution of a vec- 
tor wavepacket as it falls into an AdS geometry. The 



dynamics of the falling wavepacket are dual to the evolu- 
tion of an i?-charge jet in the same theory that we study, 
TV = 4 SYM at strong coupling. In contrast, we study 
jets of fundamental matter with large energies that scale 
like V\. 

An alternative approach to studying jets is to consider 
observables which are completely self contained and do 
not refer to a particular initial state. The best theo- 
retically motivated observables are the energy correla- 
tors introduced in [TOl HU [12] • In a state obtained from 
acting with a local current on the vacuum, the energy 
density one-point function exhibits an "antenna pattern" 
where the number of particles heading into a particular 
direction is correlated with the polarization of the cur- 
rent, giving rise to an intrinsically jet-like event. These 
antenna patterns have recently been studied in |13j for 
strongly coupled A/" = 4 (in fact for any CFT) and it 
was found that in a theory with a supergravity dual they 
are perfectly spherical and show no sign of jet-like struc- 
ture.^ A different question is to ask about the likeli- 
hood that an initial off shell photon would produce a 
quark/ antiquark pair with a particular preferred axis. 
The result of [T31 shows that an off shell photon would 
produce quark/ antiquark pairs with no preferred axis and 
on top of this no correlation between the direction of the 
quark and the antiquark. In contrast, we set up an initial 
state that has a preferred direction for both. We mostly 
focus on back-to-back jets, but as we show in Section [Tvl 
our methods apply for quarks and antiquarks heading in 
completely different directions. The results of ^13, then 
simply implies that the initial states we consider are not 
created by a local J\f = A current operator. Again, we em- 
phasize that this is not a concern for the issue of studing 
showering as the dynamics that give rise to the creation 
of a qq state via the decay of an energetic off-shell photon 
are qualitatively different in a strongly-coupled conformal 
theory than QCD, which is asymptotically free. 

The framework that allows us to do these calculations 
is the AdS/CFT correspondence |Tl[T7l[TH]. In the limit 
iVc — > oo and A > 1, TV = 4 SYM has a dual descrip- 
tion in terms of type IIB supergravity on an AdSs or 
AdSs black hole geometry times an internal 5-sphere. In 
addition, states with an energy scaling as a/A can be de- 
scribed by a classical Nambu-Goto string moving in this 
spacetime, as is implicit in the calculation of Wilson line 
expectation values of [HI HD] and was argued in detail 
in j21j . Similarly, we consider a state that consists of an 
open string which is created at a point in space at time 
t = 0. Such initial conditions locally conserve charge in 
the bulk and baryon number on the boundary. As shown 
in Fig. IT] as time evolves the string evolves from a point 
into an extended object and falls under the influence of 
gravity. The string's endpoints are charged under a C/(l) 
gauge field Am- The simplest realization of such a gauge 
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FIG. 1: A plot of a falling string at zero temperature. The 
string is created at radial coordinate m = wo and expands as it 
falls. The endpoint trajectories, shown in the figure as solid 
green and orange lines, asymptotically approach light-light 
geodesies which are shown in the figure as dotted blue lines. 



field is tiie woiidvolume gauge field of a D7 flavor brane, 
which is spacetime filling from the 5d point of view and 
wraps a 3-sphere in the internal S^. According to [H] 
the addition of such a flavor brane from the field theory 
point of view corresponds to the addition of an A/" = 2 
supersymmetric hypermultiplet in the fundamental rep- 
resentation of the SU{Nc) gauge group. 




FIG. 2: A cartoon of the bulk to boundary problem at finite 
temperature. The endpoints of strings are charged under a 
U{1) gauge field Am- The boundary of the geometry, located 
at radial coordinate u = 0, behaves like a perfect conduc- 
tor and consequently, the presence of the string endpoints 
induces a mirror current density j^ on the boundary. Via 
the AdS/CFT correspondence, the induced current density 
has the interpretation of the baryon current of a jet. Simi- 
larly, the presence of the string induces a perturbation /ia/jv 
in the metric of the bulk geometry. The behavior of the met- 
ric perturbation near the boundary encodes the information 
contained in the perturbation in the SYM stress tensor due 
to the presence of the jet. 



In the large Nc limit the coupling e^ between the gauge 
field Am and the string scales as e^ ~ ^/Nc and conse- 
quently the backreaction of the gauge field on the string's 
motion can be neglected. Similarly, in the large Nc limit 
the 5d gravitational constant k| ~ l/-^c ^o the backre- 
action of gravity on the string can be neglected too. In 
short, in the large N^ limit the string falls freely under the 
influence of gravity and the behavior of the gauge fleld 
Am (and the perturbation to the metric Hmn caused by 
the string) can be computed perturbatively. The dual 
field theory lives on the boundary of the asymptotically 
AdS geometry. The boundary behaves like an ideal elec- 
tromagnetic conductor, so the presence of the gauge field 
Am induces a current j^{x) which resides purely on the 
boundary. This is depicted schematically in Fig. [2] Via 
the AdS/CFT correspondence, this induced current has 
the physical interpretation of the baryon current density 
of a quark- antiquark pair. Similarly, the boundary be- 
havior of the metric perturbation induces a perturbation 
in the boundary stress-energy tensor T^'^. The perturba- 
tion in T'^" contains both the stress-energy of the quarks 
and the stress-energy of any emitted radiation. 

At zero temperature we show that the string endpoint 
trajectories generically approach hght-like geodesies at 
late times. While the baryon density at any point on 
the boundary may be sensitive to the behaviour of the 
source before it relaxes to geodesic motion, we are able 
to find boundary observables that are insensitive to early 
time dynamics. In particular, we find that the the center 
of charge and the angular distribution of baryon number 
are uniquely determined in terms of a single parameter 
Vs, which is the asymptotic spatial velocity of the string 
endpoint. The angular distribution of baryon number 
falls off as 7T — T7^ — sw with the angle 9 measured relative 
to the main axis of the jet. 

At finite temperature we consider string configurations 
where the string endpoints travel very far in the spatial 
directions. To generate such configurations, the strings 
must start very close to the boundary of the AdS-BH 
geometry and must have a very high energy and momen- 
tum. For such string configurations, we find that the 
endpoint motion is approximately geodesic even at early 
times. The corresponding baryon density can be made 
highly localized for a period of time and hence can be 
considered to be the baryon density of a quasi particle. 
We track the evolution of the baryon density from the 
initial event where a quark-antiquark pair was created 
to thermalization of the resulting jets and the onset of 
hydrodynamics . 

An outline of our paper is as follows. In Section K 
we write down the conventions we use. In Section III 
we discuss the dynamics of falling strings at both zero 
temperature and finite temperature. In both cases ap- 
proximate analytic solutions to the string equations of 
motion are constructed and in both cases we perform 
numerical studies of the string solutions. In Section [lll| 
we study the 5d electromagnetic problem and compute 
the baryon density induced in the boundary field theory. 



In Section IV we discuss the implications of our results 



and in Section^ we conclude. 



in the Minkowski spatial direction x = x\\. The point- like 
initial conditions are 



i(0,cr) = 0, a;(0,cr) = 0, u{0,a)=uo. 



(2.4) 



A. Conventions 

Five dimensional AdS coordinates will be denoted by 
Xm while four dimensional Minkowski coordinates will 
be denoted by x^ . Upper case Latin indices M,N,P,... 
run over 5d AdS coordinates, while Greek indices run 
over the 4d Minkowski space coordinates. Worldsheet 
coordinates will be denoted as cr° with a — 0,1. The time 
like world sheet coordinate is r = tj" while the spatial 
coordinate is a = a^. The endpoints of the string are 
located at ct = ct*. We choose coordinates such that the 
metric of the AdS-Schwarzschild (AdS-BH) geometry is 



ds' 



1/ 



-f{u) df + dx^ 



du 



2 1 



/(") 



(1.1) 



where /(u) = 1 — {u/uhY ^^id L is the AdS curvature 
radius. The coordinate u is an inverse radial coordinate; 
the boundary of the AdS-BH spacetime is at m = and 
the event horizon is located at w = u^, with T = {ttuii)~^ 
the temperature of the SYM plasma. 



The remaining necessary initial data is that of t, x and 
M at T = 0. However one of these three functions may be 
eliminated via gauge fixing. For example one may choose 
to work in the gauge t — t. We leave the detailed discus- 
sion of our particular choice of the remaining initial data 
for Section |IIB[ We simply note here that we can choose 
initial conditions such that (i) at zero temperature the 
string endpoints move infinitely far apart as i — > oo and, 
(ii) at finite temperature the initial conditions can be 
chosen such that the endpoints travel an arbitrarily large 
(but finite) distance in the Minkowski spatial directions. 
As discussed in Section HI at zero temperature strings 
whose endpoints separate to spatial infinity correspond 
to jets in the dual field theory. At finite temperature, 
strings whose endpoints can be made to travel arbitrar- 
ily far in the Minkowski spatial directions correspond to 
quasi-particle states in the dual field theory. 



A. Approximate Solutions 



II. FALLING STRINGS 

The dynamics of a classical string are governed by the 
Nambu-Goto action 



Sng — —To 



(2.1) 



where Tq = yfXjlTil? is the string tension, a and r are 
worldsheet coordinates and 7 = det 7ab with ^^h the in- 
duced worldsheet metric. The string profile is determined 
by a set of embedding functions X^^{T,a\ In terms of 
these functions we have 



and 



lah = daX ■ dbX, 



7 = (X • X'f - X^X'^ 



(2.2) 



(2.3) 



where X^^ = drX"^ and X'*^ = d„X^^ . 

The equations of motion for the embedding functions 
and the open string boundary conditions follow from set- 
ting the variation of the Nambu-Goto action to vanish. 



We write the equations of motion below in Sections II A 1 



and |II A2| The open string boundary conditions require 
the string endpoints to move at the local speed of hght 
and that their motion is transverse to the string. 

We consider initial conditions such that at worldsheet 
time T = 0, the string is mapped into a single point in 
space. We furthermore consider string profiles which lie 



At both zero temperature and finite temperature it is 
possible to obtain approximate solutions to the string 
equations of motion. Our zero temperature solution is 
valid in the t — > 00 limit. Such a limit is meaningful to 
consider at zero temperature as the string solutions will 
expand indefinitely and will correspond to jets in the dual 
field theory which separate forever. At finite tempera- 
ture, the i — > cx) limit of string solutions is not a useful 
limit to consider. This is because at finite temperature 
strings will fall asymptotically close to the black hole as 
time progresses. This corresponds to the thermalization 
of jets in the dual field theory. At finite temperature we 
perform an expansion which in essence corresponds to 
studying an arbitrarily high energy string. 

Both zero temperature and finite temperature approx- 
imations rely on the fact that there exist exact solutions 
to the string equations of motion where the string is ev- 
erywhere moving at the speed of light. Such solutions 
have infinite energy and hence are unphysical by them- 
selves, but nevertheless serve as good solutions to base 
an asymptotic expansion around. In our numeric studies 
shown below in Section |IIB| we demonstrate excellent 
agreement between our asymptotic solutions and the ex- 
act numerical solutions for a wide class of initial condi- 
tions. 



1. Zero Temperature 

To obtain the T = late time asymptotic solutions, 
we choose worldsheet coordinates such that t = t and 



a = Lp where ^ is an angular coordinate defined by 



X 

u 



£{t, if) cos if 
(.{t, ip) sinip. 



(2.5) 
(2.6) 



The domain in the {t, ip) plane in which i{t, ip) is defined 
is bounded by two curves ipsit), s = 1,2. These curves 
correspond to the motion of the string endpoints and 
are determined by imposing the open string boundary 
conditions. The open string boundary conditions require 
the endpoints to move at the speed of light and for the 
endpoint velocity to be transverse to the string. This 
implies that along the curves ipsit) 



Gm/-^^-^ = (2.8) 

at dip 

The equation of motion for i{t, ip) is 

f{f + I'^yi -f{i- i^)e" + 2e^£' cot ^ + 2i'^ cot ip 

-2fk'{£cotip + l')~£^{l-i^) =0. 

(2.9) 

Generic open string solutions to the above equation will 
fall forever into the bulk of the AdS geometry and will 
approach the speed of light. In particular, i will gener- 
ically scale like t at late times. Because we are inter- 
ested in string solutions which correspond to jets in the 
dual field theory, we seek solutions whose endpoint sep- 
aration in the spatial directions also grows like t in the 
large t limit. Such string solutions will expand at a rate 
approaching the speed of hght, and correspondingly, as 
time progresses neighboring points on string will become 
causally disconnected. Just as in inflation, any short 
wavelength structure in the string will expand and be- 
come long wavelength at late times and correspondingly, 
i'/e ^ as i ^ oo. 

Our starting point for the T = late time asymptotics 
is to expand i(t, ip) in an asymptotic expansion about 
t = oo. We write 



^(<^,i) = x('^)t + X]4(¥')i^ 



(2.10) 



Using the assumption that ^'/£ ^ as i -^ oo, one can 
conclude that x{ip) — constant. However for complete- 
ness we choose to demonstrate that this is the unique 
solution to the equations of motion which satisfies the 
open string boundary conditions. 



Consistency of the above asymptotic expansion (2.10) 



with the boundary conditions in Eqs. (2.7l-(2.8l implies 



that the functions ips(t) must have the asymptotic ex- 
pansions 



oo 
n=0 



(2.11) 



Note that this expansion does not contain any positive 
powers of t. Any positive powers of t in the expansion 
of ips{t) would imply superluminal motion for the end- 
points. 



Substituting the expansion (2.101 into the equation of 



motion (2.9) and setting the leading coefficient in the 



large t expansion to vanish, we find the following equation 
for x(<p) 

-X(l - X')x" + 2(x')' cot ip - 2xHxT (2.12) 

-(l-x')(l-2xVcot^)=0. 

Similarly, substituting the expansions into the boundary 
conditions (2.7) and (2.8) and setting the leading coef- 



ficient in the large t expansion to vanish, we find the 
following equations 



Xi^o) = 1, 
X'(^g) = 0. 



(2.13) 
(2.14) 



As advertised, the solution to Eq. (2.12) which satisfies 



the above boundary conditions is simply 

xiv) = 1- 



(2.15) 



We note that the leading order solution £ = t represents 
a null string with 7 = 0. This solution simply represents 
an arc which uniformly expands at the speed of light. 

At next to leading order in the 1/i expansion, we find 
the following equation of motion 

4^1 + £2 + 2£'^£x cot (y9 -h (^o)3 cot ^ + £!q£!^ - ^1^0 = 0- 

(2.16) 

The solution to the above equation simply relates £1 to 
£q and £\ via an algebraic equation. This feature also 
occurs at higher orders in the 1/i expansion — all €„ with 
n > 1 are determined by £q and £\ and their derivatives. 



At leading order, the boundary conditions (2.7) and (2.8) 
yield the equations 






(2.17) 
(2.18) 



The fact that the asymptotic solutions contain two es- 
sentially arbitrary (up to boundary conditions) functions 
£a{ip) and ^i(y) might seem pecuUar at first glance. How- 
ever we remind the reader that the initial conditions for a 
point-like string configuration contain two function which 
are also arbitrary up to boundary conditions. Evidently 
these two initial functions must map onto ^o((/3) and £\(^ip) 
via the equation of motion (2.9). Unfortunately we do 



not currently see any nice analytic way of relating the 
initial data for a point-like string to the functions ^o(v) 
and ^1(1^3). 

From the above discussion we see that at late times 
the strings endpoints follow the trajectories 



Xs{t) 
Us(t) 






K't, 



(2.19) 
(2.20) 



where Vs = cost^g. These trajectories are simply 
geodesies in the zero temperature AdS geometry. It is 
very easy to understand the origin of the asymptotic end- 
point geodesic behavior. After a long period of time the 
presence of the gravity implies that the radial coordinate 
of the endpoint must increase like Ug ~ t. However, by as- 
sumption are considering strings whose endpoint motion 
also scales like x^ ~ t at late times. These two conditions 
together with the condition that the endpoints move at 
the speed of light then fixes the asymptotic endpoint mo- 
tion to be a Hght-likc geodesic. 



2. Finite Temperature 

Our approach to obtaining finite temperature approx- 
imate solutions to the string equations of motion is less 
systematic than our approach at zero temperature. In 
particular we do not perform a late time expansion. 
The reason is simple — at asymptotically late times the 
string endpoints always fall into the black hole, and con- 
sequently the endpoint motion in the spatial directions 
ceases. This corresponds to the energy loss and thermal- 
ization of jets in the dual field theory. The late time be- 
havior of any string solution at finite temperature simply 
consists of a string of constant spatial extent, whose ra- 
dial coordinate asymptotically approaches the event hori- 
zon. 

Instead of performing a i — > oo expansion, we construct 
an asymptotic expansion for a string solution with very 
high energy. The validity of our expansion will turn out 
to be equivalent to the statement that the energy density 
and energy flux on the string are parametrically high. 
From the field theory point of view this is not surprising. 
The quark that the string is describing loses energy to the 
SYM plasma at a rate dEq^ai^/dx ^ E" for some n > 0. 
Therefore asymptotically high energy states will deposit 
a large amount of energy per unit length into the plasma. 
From the gravitational point of view this corresponds to 
having a large energy flux and energy density along the 
string. 

Generically the endpoint motion of a falling string is 
bounded by that of a light-like geodesic which starts at 
u — uq. It is easy to work out the equations of motion 
for light-like geodesies in the AdS-BH geometry. They 
read 



we see 



dxa 



dt 

du, 



geo 



dt 



Pie 



/) 



e 



(2.21) 
(2.22) 



where ^ is a constant which specifies the geodesic. Note 
that for any geodesic which obtains a minimal radial co- 
ordinate Uniin on its trajcctory, the para meter ^ must 
satisfy ^ > /(umin)- Dividing Eq. (|2.2l| by Eq. (|2.22[) 



dx„ 



du 



e-f 



(2.23) 



Eq. ( |2.23 l may be integrated to determine the total 
spatial extent that the geodesic travels. In the limit 
e -^ f{uo) and uq — > 0, this distance scales like u\/uq. 

We consider string configurations where the spatial 
component of the string endpoint motion Xg p^ t for times 
t <C u\/uo. Under this condition the endpoint velocity 
in the Minkowski spatial directions will be very close to 
the local speed of light. Because open string endpoints 
must always travel at the speed of light, the velocity in 
the radial direction must therefore be small and corre- 
spondingly, the radial coordinate of the string endpoints 
must change very slowly for times t <^ u\/uq. 

Just as at zero temperature, the assumption that the 
Xg K t implies the string must stretch and expand as 
time evolves. For sufficiently smooth initial conditions^ 
this implies short wavelength perturbations in the ini- 
tial structure of the string will expand and become long 
wavelength, resulting in a smooth string profile at late 
times. Moreover, as the string endpoints separate, the 
middle of the string must fall rapidly toward the event 
horizon. In our numerical studies discussed in Section 



IIB| we have found that this typically occurs over a time 
scale /S.t ^ Uh- This is the time scale for a geodesic with 
f <C 1 to fall into the horizon. We see this for both 
of the numerically generated string solutions in Figs. [6] 
and [9J The origin of this behavior can be understood 
as follows. Consider the string at times shortly after the 
creation event. It will have expanded to a size r^ t. By 
construction the string will have a very large momentum 
density in the spatial directions. The momentum den- 
sity must be very inhomogeneous so that each endpoint 
can move off in its own direction. As time progresses 
the parts of the string with the highest momentum den- 
sity will evolve to be located near the string endpoints. 
Portions of the string with low momentum density will 
lag behind the endpoints and fall relatively unimpeded 
toward the horizon. After the middle of the string falls 
into the event horizon, the trajectories of the two halves 
of the string which extend up from the horizon will be 
uncorrelated. As time progresses, each half of the string 
will translate along in the spatial directions and the end- 
point will slowly fall. We see this nicely for the numerical 
strings in Figs. |6]and[9] 



^ Initial conditions where there is an arbitrarily large amount of 
structure in the string at early times do not have to relax to a 
smooth string profile after any given amount of time At. While 
the initial structure will be infiated as time progresses, because 
the string endpoint can only travel a distance ~ u^/uq, one can 
always cook up initial conditions such that the initial structure 
never relaxes in this time interval. We do not consider such initial 
conditions in this paper. 



With the above motivation in mind, our starting point 
for the construction of approximate finite temperature 
string solutions is to consider a half string which asymp- 
totically approaches the event horizon at spatial dis- 
tances far from the string endpoint. We look for solutions 
which approximately translate at constant velocity in the 
spatial directions. These solutions will be approximately 
stationary in the rest frame of the string endpoint during 
the time interval uu-^t <^ "^/"o- 

We choose worldsheet coordinates t = t, a — u. In 
these coordinates the embedding functions are deter- 
mined by one function x{t,u). The domain in the (t,u) 
plane in which x{t, u) is defined is bounded by a curve 
Us{t). This curve corresponds to the trajectory of the 
string endpoint and is determined by imposing the open 
string boundary conditions. The equation of motion of 
x{t, u) is 

2u (1 + fx'^) X - 2uf (/ - i^) x" - iufx'xx' (2.24) 
-I- 4/ (2 - / (1 - x'^)) x' - 4 (3 - 2/) x'x^ = 0. 



We seek a solution to Eq. (2.241 of the form 

x{t,u) = Xstce,dyit,u) + exi{t,u) + 0{e^) (2.25) 



where 



a;stoady (t, U) = C^ + Xq (u) , 



(2.26) 



is a steady state solution to the equations of motion with 
^ a constant and e is a formal expansion parameter. Sim- 
ilarly, we write 



4t)^u^o{t) + eulit)+Oie^). 



(2.27) 



The function xi{t,u) characterizes the perturbations in 
the string which have inflated to long wavelengths. 

We note that having nontrivial time dependence in 
Us{t) does not contradict the steady state assumption in 
Eq. (2.261. As we shall soon see, the steady state condi- 



tion implies that uo{t) must satisfy a geodesic equation. 
At leading order in e the equation of motion (2.241 lead 
to the following equation of motion for xo{u) 



+A[{2-f)f-e{i-v)\^ 



(2.28) 



0. 



The general solution to the above equation is given by 
functions which satisfy 



du ) ulp (1 - CJ) 



(2.29) 



where C is an arbitrary constant. Before integrating 
Eq. (2.29), we first determine C by imposing the open 
string boundary conditions. The reader may note that 



Eq. ( 2.29[ ) is the equation of motion for the dragging 
string of Refs [HIIT]- There, the string obeys a Dirichlet 
condition at the endpoint and consequently C is found 



to be 1/C^. We will find a different value of C because 
we solve the problem for different boundary conditions. 
As a result, we generally obtain a different solution than 
the string of Refs [6,7. 

With our choice of worldsheet coordinates, the open 
string boundary conditions are simply 

Gmn'^-^'^-^ = 0, (2.30) 



G 



dt dt 



MN- 



0, 



(2.31) 



dt du 

where all quantities are evaluated along the curve Ug(t). 
At leading order in e these boundary conditions lead to 
the two equations 

2 

(2.32) 



dxQ 
du 

duf) 



e-f 
Pie- 



f) 



e 



(2.33) 



where all quantities are evaluated along the curve UQ(t). 
Comparing Eq. (2.321 with Eq. (2.291, we see the two 
equations agree if C ~ 1. Furthermore, comparing 



Eq. (2.33) with Eq. (2.22), we see that to leading order 



in e, the string endpoints follow light-like geodesies. 
With C = 1 Eq. (|2.29[) may be integrated. Taking the 



negative root of the equation so that the string profile 
trails the endpoint, we can find the profile in terms of 
Appell hypergeometric functions. For the special case of 
f ^ 1, xo{u) reduces to 



xo{u) 



Uh 

2 



_i U 1 Uh ~ u 

tan 1 — log 

Uh 2 Uh + u 



(2.34) 



which is the well know trailing string profile of Refs |6l [7] . 

Just as at zero temperature, the leading order solu- 
tion to the string equations of motion is a null string 
with 7(xstoady) = 0. Including the 0{e) correction to the 
string profile implies ^{x{t,u)) = 0{e). The energy den- 
sity and energy fiux of the string both are proportional 
to l/\/— 7. Therefore the validity of our asymptotic ex- 
pansion implies the energy density and energy fiux of the 
string are very high. 

It is straightforward to work out the equations of mo- 
tion for the perturbation xi(t,u). We do not carry out 
this exercise here but we simply note that one can show 
that the perturbations do not become large as time pro- 
gresses. In particular if e is small, then exi{t,u) remains 
small for all times. We therefore expect the perturbative 
expansion in Eq. (2.25 1 to be valid outside the time win- 
dow Uh ^ t <^ ujJuq. Indeed as we discuss in Section 
|IIB| the agreement between the steady state result and 
our numerical results agrees for all times i ~ a few Uh- 



B. Numerical String Solutions 

In addition to studying the asymptotic behaviour of 
falling strings via the expansions in Section II A we can 



numerically solve the string equations of motion to find 
the worldsheet for arbitrary initial conditions. These so- 
lutions allow us to test the accuracy of the expansions in 
Section [II A| and to study string behaviour in the regimes 
where the expansions do not hold. 

For numerical reasons discussed in Ref [6], we have 
found it convenient to work with the Polyakov action 
instead of the Nambu-Goto action. The Nambu-Goto 
action is classically equivalent to the Polyakov action 



Sv 



' 2 



drrdax 



-r,r,'^''daX^'dbX''GMN, 



(2.35) 



where r]ab is the worldsheet metric. The equations of 
motion follow by varying of the action with respect to 
X^' and r]ab- The variation with respect to X^^ yields 
the equation of motion 



da [V^V^'CMNdbX^] = 



-daX^'dbXP, 

(2.36) 



2 V •'•' QXM 

together with the open string boundary condition 

<,(T,a*) = 0, (2.37) 

where 

(55p 



</(T,cr) = 



= -nV^ri'^'GMNdaX 



N 



5X'M{T,a) 

(2.38) 
and a = a* at the string endpoints. 

The variation of the action with respect to rjab yields 
the worldsheet constraint 



1 



lab 



VabV led- 



The constraint equation implies 



(2.39) 



(2.40) 



which when substituted back into the Polyakov action, 
easily shows the classical equivalence to the Nambu-Goto 
action. 

One may fix the worldsheet coordinates (r, a) by mak- 
ing a choice for the worldsheet metric rjab- Following Ref 
[B], we choose to fix 



iVab) 



— E(a;,w) 
1 



S(a:,M) 



(2.41) 



where S(2^, u) is a stretching function, which we take to 
be a function of x{t, a) and u(t, a) only. With this choice 
of worldsheet metric, the worldsheet constraints (2.39) 
read 



X^- 



X-X' ^0, 
Y^X'"^ = 0. 



(2.42) 
(2.43) 



Since we choose to study point-like initial conditions, we 
must choose initial data X*^(0,cr) consistent with the 



cons traint equation Eq. (2.431 and boundary conditions 
Eq. (2.371. This is generically satisfied by specifying t in 
terms of x and u via 



ft' =x' + 



f 



in addition to requiring the boundary conditions 
i'{a*)^x'{a*)^u{a*)=Q 



(2.44) 



(2.45) 



at r = 0, where a* denotes the endpoints of the string. 

To choose the functions i(0, a) and u(0, a) as well as 
the initial radial coordinate uq, we must consider what 
kind of states we wish to study in the field theory. We 
first consider the zero temperature theory. As we discuss 



below in Section III the endpoints of the strings roughly 
correspond to the locations of the jets in the boundary 
field theory. Field theory states with jets that propagate 
to infinity therefore correspond to falling strings whose 
endpoints escape to infinity in the spatial directions. We 
therefore choose to study string initial conditions with 
enough energy and momentum so that the endpoints sep- 
arate from each other for all times. In the finite temper- 
ature theory, the distance the endpoints can separate is 
bounded. The field theory states here have a natural 
quasiparticle description until the jets thermalize. We 
therefore choose to study initial conditions where we can 
tune how far the jets travel, i.e. how long the quasipar- 
ticle lives before thermalizing. 

Beyond these restrictions, the initial conditions are not 
further constrained by physical requirements. Different 
initial conditions correspond to different structures of the 
initial string. What we want to be able to do is to iden- 
tify aspects of the resulting jets that are universal, but 
also aspects that depend crucially on the choice of initial 
conditions. 



1. Summary of Results 

We have examined a number of initial conditions in 
the zero and finite temperature cases and in each case 
have confirmed the approximations in Section |II A[ In 
particular, the string endpoints asymptote to geodesies 
in each geometry and fiuctuations along the string at 
early times are infiated in wavelength as the string falls. 

We will now showcase four representative solutions — 
two at zero temperature and two at finite temperature. 
The details of how these strings are generated are given 
below in Section |IIB 2 The first string we examine is 
plotted in Fig. |3] with the null string solution I — t over- 
laid. From the figure, we see that the the difference 
between the null solution and the numeric solution be- 
comes small compared to the size of the string as time 
progresses. Moreover, we see that the endpoint trajec- 
tory is almost completely geodesic. 

The above string was generated with rather symmet- 
ric initial conditions. To better test the series solution 
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FIG. 3: A falling string at zero temperature. The solid 
blue curve represents the string at three successive times 
ii = 15uo, t2 = 30iio, and fa = 45uo. The dashed red curve 
represents the null string approximation £ — t a.t the same 
three times. The orange and green solid curves represent the 
trajectories of the string endpoints while the dotted blue lines 
represent the geodesic fit the the endpoint trajectories. The 
num erical string was generated with the initial conditions in 
Eq. (2.481 at the large blue dot near the top of the plot. As 



time progresses, the difference between the null string and the 
numeric string becomes small compared to the overall size of 
the string. 



FIG. 4: An asymmetric string at zero temperature generated 
from the initial conditions Eq. (2.481. The solid blue curve 
represents the string at three successive times fi — 7uo, i2 = 
14wo, and is — 21uo. The dashed red curve represents the 
null string approximation i = t a.t the same three times. The 
orange and green solid curves represent the trajectories of 
the string endpoints while the dotted blue lines represent the 
geodesic fit to the endpoint trajectories. As time progresses, 
the perturbations in the string profile inflate and become long 
wavelength. Correspondingly, the difi'erence between the null 
string and the numeric string becomes small compared to the 
overall size of the string. 



Eq. (2.101, we therefore consider a second zero temper- 
ature string with large fluctuations at early times. Such 
a string is shown in Fig. |4] The profile has two large 
fluctuations at early times that grow in size as the string 
falls. At these early times, the circular arc solution is not 
a good flt to the string profile; there is a large perturba- 
tion. However as time progresses the short wavelength 
structure in the perturbation inflates to long wavelengths 
and the size of the perturbation becomes small compared 
to the growing size of the string. Again from the fig- 
ure, one can see that the endpoints approach light-like 
geodesic trajectories. 



According to the expansion given in Eq. (2.101, the dif 



ference between any expanding string and the null string 
i = t should have the late time form 



i{t,^)-t = eo{^) + 



hiv) 



(2.46) 



so for fixed (p, {(.{t, ip) — t)t should be linear in t as i — j oo. 
This observation provides a simple test of the form of 
the asymptotic expansion given in Eq. (2.101. We plot 
{£{t, if) — t)t in Fig. Is] for the above asymmetric string 
configuration shown in Fig. [4] As is evident from the 
figure, for fixed tp, {£{t, ip) — t)t is linear at late times, 
reinforcing the form of the expansion given in Eq. (2.10 1. 



In the same spirit as with the zero temperature plots 
above, we exhibit two finite temperature strings — one 
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FIG. 5: A plot of the deviation {l{t,ip) — t) f of the string 
shown in Fig. [4] from the null string. According to the asymp- 
totic expansion given in Eq. (2.101, for constant values of ip, 
{£{t, (fi) ~ t)t should be a linear function of t at late times. 
The linearity of the plots at late times reinforces the form of 
the expansion in Eq. (2.101. 



smooth string solution and one solution with large fiuc- 
tuations at early times. A smooth string is shown in 
Fig. [6J It is a symmetric state about x = whose end- 
points travel a spatial distance of 19.25/7rT. The dotted 
blue line in the plot represents a geodesic fit to match the 
stopping distance of the string endpoint. The dashed red 
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curve shows the corresponding null string plotted on top 
of the numerical string. As is evident from the figure, the 
null string fits the numerical solution very well for times 
i ^ a few Ufi- 

To see how well the endpoint motion is approximated 
by a geodesic, we plot f{dxs/dt)^^ in Fig. m For a 
geodesic, this is the constant ^. The geodesic fit to the 
stopping distance is plotted in addition to the endpoint 
data. In the plot, we see that f{dxs/dt)~^ changes by 
roughly one part in ten thousand over the course of the 
trajectory of the endpoint. This is completely consistent 
with the expansion given in Eq. ( |2.27 l. In particular, 
the value of f{dxs/dt)^^ should change by an amount 
parameterized by e. Evidently, in this case e ~ 10^**. 

We also exhibit an asymmetric finite temperature 
string, as shown at several times in Fig. [9l A close-up 
of the early-time profile is plotted in Fig. |8] As with the 
asymmetric zero temperature string in Fig. |4] the early- 
time fluctuations quickly inflate to long wavelengths. The 
string quickly relaxes over a time t ~ u/j to match a null 
string with fluctuations, matching the predictions of Sec- 
tion |llX2 



2. Numerical Details 

We generate our numerical strings by solving the string 
equations of motion with a canned PDE solver: Mathe- 
matica's NDSolve. We do so with the same technology 
used in Ref [6J . The crucial difference with the numerical 
work in Ref [5] is that we must optimize the stretch- 
ing function on a case-by-case basis. We essentially tune 
the worldsheet metric to each solution so that gradients 
along the string are minimized at all times in the string's 
evolution. We do so with a stretching function of the 
form 



E{x,u)^{x^ + iy 



1 



1 -uo 



fu^Y 



V u 



(2.47) 



Our zero temperature strings can generally be tracked 
with high accuracy if we only turn on the ^ part of 
this stretching function. Both of the zero temperature 
strings plotted in this paper simply had p — 2 and in each 
case we were able to track the string down to distances 
past IO^wq- The configuration exhibited in Fig. ^ was 
generated with the initial conditions 

X — 4uocoscrexp [— i(l — cos^ cr)] , (2.48a) 
ii = 2uo, (2.48b) 

where we have specified the endpoint locations a* to be 
and TT for all of our numerical strings. 

Here, we would like to give the reader some intuition 
about why we chose these initial conditions by highlight- 
ing their salient features. First, the exponential factor in 
X allows us to control how much momentum and energy 
is initially located near the string endpoints. This should 
roughly correspond to how mvich momentum and energy 



are stored near the quark and antiquark at early times. 
For this solution, the string therefore has more (but not 
much more) energy near its endpoints than in the mid- 
dle. The rest of the initial conditions are chosen to be as 
simple as possible such that the endpoints fiy apart from 
each other at t = 0. 

The asymmetric string shown in Fig. [4] was generated 
with a more complicated set of initial conditions given 

by 



-4?io (cos a + h cos 3(j) , 



u = Uq 



2 cos 4(7 — cosSct) 



(2.49a) 
(2.49b) 



Here we chose initial conditions so that the initial veloci- 
ties along the string, dx/dt — x/t and du/dt — ii/t, would 
have large fluctuations at early times. These evolve to the 
bumps that inflate in size as the string falls. 

Similar to the first zero temperature solution, the 
initial conditions for the symmetric finite temperature 
string in Fig. |6]were 



Mo 

X 



h^h, 



626^0 COS cr exp [—16(1 
u = uoV/("o)(l -cos2cr). 



(2.50a) 

cos^ a)] , (2.50b) 

(2.50c) 



In this case, we needed to turn on the other parts of the 
stretching function Eq. (2.471. For this string, we found 
the choice m — 0.6, n — 1, and p = 2 to be effective. 
Also, the initial conditions here are somewhat more ex- 
treme than in the zero temperature cases above. Our 
choice reflects an attempt to make a long-lived string by 
(i) creating the string near the boundary, (ii) putting 
most of the energy and momentum near the endpoints, 
and (iii) giving the endpoints zero initial radial veloc- 
ity. This last ch oice al lows the string endpoints to travel 
along an ^ « \//('Uo) geodesic for a significant fraction 
of their trajectory, whereby they travel much further in 
the spatial direction than an ^ = 1 geodesic. See Section 



IV for more details on this point. 




The initial conditions for the asymmetric 


finite tem- 


perature string in Fig. [9] were 




Wo = jqUu, 


(2.51a) 


X = 30uo (cos cr -1- ^ cos 3(7 - ^) , 


(2.51b) 


ii — uo\/f{uo) (5 — 2 cos 4(7 — |cos5cr) 


. (2.51c) 



We used a slightly different stretching function for this 
string with m = 0.7, n = 1.5, and p = 2. Here, the initial 
conditions are similar to those used for the asymmetric 
zero temperature string of Fig. |4] in that we generate a 
string that has large fluctuations in the velocities dx/dt 
and du/dt at small times. We also create the string near 
the boundary with a large amount of energy near the 
endpoints so that, like the long-lived finite temperature 
string in Fig. [91 the endpoints will travel some long dis- 
tance before falling into the horizon. 



11 



f 



0.5 



1 



1 



5 10 



15 



FIG. 6: A symmetric string at finite temperature generated from the initiai conditions Eq. (2.501. The string, shown as the 
soiid blue curve, is shown at successive times t\ = l/5nT, fa = 2.3/nT, ta = S.S/ttT, and t4 = IS.S/ttT with the corresponding 
null string, shown as the dashed red curve, plotted at the same times. The solid green and orange curves represent the endpoint 
trajectories while the dotted blue curve represents the geodesic fit to the endpoint trajectory. The null string agrees very well 
with the numeric string at times i ~ a few Uh- The geodesic fit to the endpoint trajectory was obtained by matching to the 
total distance traveled by the endpoint. 
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FIG. 7: The quantity f{dxs/dt)^^ for the string of Fig. [mis plotted with the solid blue line as a function of time. For a 
geodesic, this is the conserved quantity ^. We also plot the geodesic fit to the endpoint trajectory with the dashed purple line: 
^ = 0.999994. As is evident, f{dxs/dt)~^ is approximately constant over the trajectory of the st ring endpoint. The deviations 
of f{dxs/dt)~^ from a constant give a measure of the expansion parameter e given in Eq. (2.271. 



III. THE ELECTROMAGNETIC PROBLEM 
AND THE BOUNDARY BARYON DENSITY 

In the large Nc limit the action for the U{1) gauge 
field living on the D7 brane reduces to that of ordinary 
curved space Maxwell equations. We therefore start with 
the Nc ^ oo effective action 



S — SeM + 'S'int ■ 

The electromagnetic action is given by 



(3.1) 



where J^ is the current corresponding to the string end- 
points. 

The equations of motion for the gauge field follow from 
taking the variation of the action (3.1 ). The variation of 

5 is 



SS = / d^xV^GSAN { ^DmF"'^^ - J^ 



(3.4) 



(3.2) 



where Fmn is the field strength tensor corresponding to 



the gauge field Am and ^ 
is given by 



47r2 



The interaction action 



S'int = / SxV^J-A. 



(3.3) 



The first term in Eq. (3.4 1 gives the equations of motion 
for the gauge field 



1 



DmF'^'' ^ J^, 



(3.5) 



while the surface term gives the boundary baryon current 



12 



"h 




FIG. 8: A closeup of the string in Fig. [9]at time t = I/SttT. 
This plot shows the initial structure in the string which is 
inflated to long wavelengths at late times. 



conditions. The differential operator in (3.9 1 has singu- 
lar points at It = and u = Uh with exponents 1 and 
±icjM;i/4. Regularity at the boundary implies g<^{u) ^ u 
as u ^ while the requirement that the black hole not 
radiate [23 implies that g>{u) ~ (m— u/i)""^"''/'* near the 
horizon. The overall normalization of (?< (u) may be fixed 
by requiring Irm^^Q g^(u)/u = 1. From Eq. (3.7 1 we 



therefore find that the boundary charge density is given 
by3 

p{lu, q)^^ r du' Giu, q, ^')^5(c^, 9, u') (3.14) 



p2 / 

e Jo 



where 



g{uj,q,u) 



gy{w,q,u) 

W{uj,q,u) 



(3.15) 



is the Fourier space bulk to boundary propagator 



density 



ss 



1™ ^ . ,r Tj (3-6) 

-^ lim v/-G'(t,a;,u)i^^^(i,a;,M). (3.7) 



e^ u^o 



We note that the above expression for the boundary cur- 
rent density can also be obtained by integrating the equa- 
tions of motion (3.51 over a gaussian pillbox enclosing the 
boundary. 



From Eq. ( 3.7 1 we see that the boundary charge density 



is determined by the near boundary behavior of the radial 
component of the electric field 



E^ = Fi 



05- 



(3.8) 



The equations of motion for _E5 can easily be worked out 
from Eq. (3.51. Introducing a spacetime Fourier trans- 



form, we find the following equations of motion for the 
mode amplitudes E^(uj,q,u) 



where 



^5 = 



^s^-.-s + B^Er, — Si 



(3.9) 



uf ' 




(3.10) 


P 


uf-f 


(3.11) 


e^L^ (2Jo - 


- iuujj^ — uJq) 


(3.12) 



We may solve Eq. (3.9 1 with a Greens function G{u, u') 
constructed out of homogenous solutions 

G{u,u') = g>{uy)g^{u^)/W{u') (3.13) 

where W(u) is the Wronskian of g^ and gy. The appro- 
priate homogenous solutions are dictated by boundary 



A. Zero Temperature Baryon Density 

At T = the AdS radial coordinate ranges over 
the interval (0,oo) and / = 1. In this case homoge- 
nous solutions to Eq. (3.9) may be obtained analyti- 
cally. The boundary condition on gy{u) at u = cx) is 
limu^oo g>{u) = 0. We find 



9>{u) = uKoiuy/q'^ -^2), 



(3.16) 
(3.17) 



where Io{z) and Ko{z) are mod ified Bessel functions. We 
note that as given in Eq. (3.171, g>{u) has an ill defined 
M — > cxD limit. This can be remedied by giving uj an 
infinitesimal imaginary part. Causality dictates that the 
retarded Greens functions be analytic in the upper half 
uj plane, so the appropriate prescription is to send lo — > 
Lu+ie. With this prescription g>(u) vanishes at u = oo for 
all q and co. The resulting bulk to boundary propagator 
is therefore 



g{uj, q, u) = -Ko{u^/q^ -uj'^). 



(3.18) 



Using the aforementioned ie prescription, Eq. (3.18) may 



be Fourier transformed to real space to obtain 

g{t,x,u) = -0{t)^Siw), (3.19) 

TT aw 

where w is defined to be 



w = -t^ + x^ + u^. 



(3.20) 



The zero temperature real space source for E^ is 



S^^e^L^[-d^{Ja/u^) + doJ5/i 



(3.21) 



We note that 55 <x ^ so the baryon density is independent of 
both e and L as it must be. 
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FIG. 9; An asymmetric string at finite temperature generated from the initial conditions Eg. |2.5l| The solution and the fitted 
null string is shown in the same manner as in Fig. [6] at times i = I/SttT, I.T/ttT, 'i/ixT, and %.! J-kT. 



Combing the real space version of Eq. (3.141 as well as 



Eq. (3.19 1 and (3.21 1 and integrating by parts, we obtain 



p(t,x) = I? du'd'^x' du'Q{t-t' ,x~x' ,u')jQ{t' ,x' ,u') 



+ doGit-t' ,x-x' ,u)J5{t' ,x' ,u) 



,/2 ■ 



(3.22) 



Using coordinate time as a parameter for the string end- 
point trajectories, the current corresponding to the string 
endpoints is 

J^(^,a;,u) = gft)^(-1)- ^^ 6^x^Xs)S{u-Us) 

(3.23) 
where the endpoint trajectory is given by X*^(t) and the 
factor of (— 1)** comes from the fact that the string end- 
points ar e opp ositely charged. Substituting Eq. ( |3.23 ) 
into Eq. (3.221 and integrating over x' and u' , we obtain 



Pit, x) = ^(-1)^ j dt'^S{w,)T,it, t'). (3.24) 



where 



r,(i,t')- — eit-t')9it')us{t') 



Us{t') + {t-t')^ 



dt' 

(3.25) 



and Ws is evaluated at 

Ws = -{t - t'f + {x- Xs{t')Y + {us{t')f. (3.26) 



1. Late Time Asymptotics 

In general the value of the baryon density at a partic- 
ular point (t, x) depends on the history of the sources in 
the bulk of the AdS geometry. However, the dependence 
on the history of sources can be reduced if one consid- 
ers quantities which are averaged in some sense, such as 



moments of the baryon density or the angular distribu- 
tion of baryon number discussed below. Moreover, the 
fact that the trajectories of the string endpoints asymp- 
totically approach light-Hke geodesies suggests that the 
late time behavior of some averaged quantities should be 
independent of the history of the the bulk sources. We 
explicitly demonstrate that this is true for for the first 
moment of the baryon density and for the angular distri- 
bution of baryon number defined below in Eq. (3.341. 



Consider the first moment of the baryon density asso- 
ciated with jet s, the baryonic center of charge position 



Xsit) 



J d^xx ps{t,x) 
J d^xps{t,x) 



(3.27) 



where ps{t,x) is the baryon density corresponding to a 
single string endpoint labeled by s. The denominator is 
simply J d'^x ps{t,x) — (—1)*+^. This follows from the 
fact that the mirror charge induced on the boundary is 
always opposite of the corresponding source charge in the 
bulk. From Eq. (3.241 we therefore have 



Xs{t) 



dt'd^xx-—rS(ws)TJt,t'). (3.28) 

dwi 



Shifting integration variables and performing a couple of 
integrations by parts, we obtain 



x,{t)^n I dt' — 



d 



x,{t')Ts{t.t') 

Zs{t,t') 



d{Zs{t,t')) 

^/zjt;t^) 



where we have defined 



{t~t'f~{uS')Y 



(3.29) 



(3.30) 



and Z, = dZs/dt'. 

In the limit t —^ oo the domain of integration in 
Eq. (|3.29| will be dominated by t' ~ t. As dis- 



cussed in Section III A II the late time behavior of the 
string endpoint motion asymptotically approaches light- 
like geodesies. The particular trajectories are given in 
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Eqs. (2.19)-(2.20). When substituted into Eq. (3.291, and 



these trajectories yield the late time behavior 

Xs{t) = Vst 

where 



1 



1- 



1 

V? 



tanh ^ Vs, 



(3.31) 



(3.32) 



and Vs is the asymptotic velocity of the string endpoint 
in the Minkowski spatial directions. (The directional de- 
pendence of Vg is the same as that of Vs.) We therefore 
see that the velocity of the center of baryon density does 
not generically equal the velocity of the the string end- 
point in the spatial directions. However in the Vs —^ I 
limit the two velocities do agree. In particular 



ws = l- 



1 -log 



l-Vs 



il-Vs) + Oiil-Vsf). (3.33) 



In this limit the string endpoint trajectory approaches a 
line of constant radial coordinate u and the position of 
the center of baryon density coincides with the spatial 
coordinate of the string endpoints. 

At zero temperature it is useful to consider the angular 
distribution of baryon number 



dBs 

dn 



r'^drps{t, x) 



(3.34) 



where r = \x\. This quantity is in spirit similar to the jet 
energy profiles used in collider physics, see e.g. [24ll25]. 
For localized lumps of baryon density which escape to 
spatial infinity, this function measures the angular dis- 
tribution of baryon number of the jet.^ 
From Eq. (|3.34|) and Eq. (|3.24|) we have 



dBs 
dfl 



i-iy 



r^drdt'^5{ws)Ts{t,t'). 
dwi 



(3.35) 



Using the delta function to carry our the r integral, we 
obtain^ 



dBs 

dn 



(-1) 



s+l 



^^, As{t,t')-3xs{tTcosH ^ 

8As{t,t')l ^^'^ 



(3.36) 
where 6 is the angle between Xg and x and 
As{t, t') = {t- t'f - (1 - cos^ e) Xs{t'f - Us{t'f (3.37) 



* Alternatively, one could measure the angular distribution of 
baryon number at infinity with the integrated baryon number 
flux r^ f dtx ■ j(t, x). 

^ Technically we should be summing over several zeros of the delta 
function in Eq. | |3.35[ l. However if one assumes cos 6 < 0, then 
there is only one zero of the delta function with r > 0. We con- 
sider cos 9 < and obtain the general result for —jjf- by analytic 
continuation. 



<i>s{t, t') = r,(t, t')e{{t ~ t'f - xsit'f - usit'f) . (3.38) 



In the t ^ oo limit the integration in Eq. ( |3.36 l will 
be dominated by times t' ^ t. Just as for the center of 
charge, we may therefore use the asymptotic endpoint 
trajecto ries in Eq. (|3.36 ). Su bstitu ting the trajectories 
in Eqs. ([2l9l)-(|2.20 into Eq. (plel, we find 



dBs ^, .s+i^-Vs' 
dn ^ ' An 



dt 



,tt'{2U'+2t'^Vs^cos''0-e) 
(3.39) 



Carrying out the integration in Eq. (3.39), we find the 
remarkably simple result 



dn 47r(l-V;cos( 



n2 ■ 



(3.40) 



From Eq. (3.401 we may compute the average opening 
angle of the baryon density. A simple calculation shows 
that for each jet the opening angle is given by 



TT 

2K 



Vi - v;2 -1 + 14 



In the K 



1 limit this becomes 
- irVT^Vs 



V2 



+ Oil-Vs). 



(3.41) 



(3.42) 



which shows that jets can be made highly coUimated and 
focused. Moreover, from Eq. (3.331, we see the limit 



where jet are focussed coincides with the limit where the 
center of baryon density velocity is close to the speed of 
light. 



B. Finite Temperature Baryon Density 

At finite temperature we obtain the boundary baryon 
density by numerically solving the Maxwell equations 
(3.9). For a given momentum (w,q), the source ^5 de- 



fined in Eq. (3.12 1, is computed by using the electromag- 



netic currents corresponding to the numerical string solu- 
tions discussed in Section II B In addition, for each {uj, q) 
the homogeneous solution 5>(u) is evaluated by numer- 
ically integrating the homogeneous differential equation 
(3.9) (without sources) outward from the horizon. Then 



the solution g<,{u) is evaluated by numerically integrat- 
ing the same equations (without sources) inward from 
the boundary. Given these numerically determined ho- 
mogeneous solutions, the baryon density is evaluated by 
numerically performing the radial integrals in Eq. (3.14), 



with the source S^ . The subsequent Fourier space baryon 
density is then numerically Fourier transformed to real 
space. In what follows, we consider the string generated 
by the initial conditions given in Eq. (2.50). 
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FIG. 10: A snapshot of the string profile and the corresponding baryon density at an instant of time well before thermalization. 
The black plane represents the event horizon of black hole and the surface on the top represents the boundary baryon density. 
The red curve is the string and the vertical axis in the plot serves as both the AdS radial coordinate and the magnitude of the 
baryon density. The baryon density is highly localized above the string endpoints. 



Fig. [To] shows a snapshot of a typical finite tempera- 
ture string/baryon density configuration at a time before 
thcrmahzation, but much after the initial event where the 
string was created.^ The black plane at the bottom of the 
figure represents the event horizon of the black hole and 
the surface on the top of the figure represents the bound- 
ary baryon density. The red curve is the string and the 
vertical axis in the plot serves as both the AdS radial 
coordinate and the magnitude of the baryon density. As 
is apparent from the figure, the baryon density is highly 
localized above the string endpoints. This type of be- 
havior is universal at times t such that u/j <C i <C u\/uq 
where uq is the initial radial coordinate where the string 
was created. 

The fact that the baryon density is localized above the 
string endpoints is easy to understand from the zero tem- 
perature analysis in the previous section. First of all, in 
the limit where the sources in the bulk are very close 
to the boundary, the electromagnetic bulk to boundary 
propagator may be approximated with its zero tempera- 
ture limit. During the time interval Uh ^ t <^ uf^/ug, 
the radial coordinate of the string endpoint will have 
not changed much from its initial value of uq. In par- 
ticular, as discussed in Section [ll A 2| the trajectories of 
the string endpoints are approximately given by light-like 



geodesies in the AdS-BH geometry, so from Eq. (2.231 



the rate that the endpoint falls toward the black hole 
will be suppressed by some power of u^/uh-'^ The source 
^5 for the radial component of the electric electric field 
Ec, will contain temperature dependence — at leading 
order in u^/uh this temperature dependence is simply 
the temperature dependence of the trajectories. There- 
fore, to leading order in uo/uh and during the window 
Uh ^ t <^ u\ /uq, we may use the zero temperature 
expression Eq. (3.311 for the center of the baryon den- 



sity but with the finite temperature trajectories for the 
string endpoints. For string configurations whose end- 
points propagate very far in the spatial directions, to 
leading order in uo/uh we have Xs(t) = t d uring the win- 
dow Uh ^ t <^ u\/uq. From Eq. (3.311 we therefore 



conclude that Xs{t) — t as well, showing that the po- 
sition of the string endpoints agrees with the center of 
baryon density. 



To elaborate on this behavior, we plot in Fig. 11 the 
dipole moment of the baryon density and the difference in 
spatial coordinates of the string endpoints. As is evident 
from the figure, the two curves agree very well over the 
vast majority of the string endpoint trajectories. How- 
ever while the figure demonstrates that as t ^ oo the 
dipole moment agrees with the endpoint separation, at 



See |http : // www . phys . Washingto n ■ edu/""karch/Papers/| 

Falling/3D.gif Ito download a movie of this event. 



"^ For example, for the special case of ^ = 1, Eq. (|2.23|l implies 

2 

that Us{t) ~ — §- in the window Uf^ <^t <^ u^/uq. 
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times t 



\/uq, there is a small discrepancy between 



dipole moment and the endpoint separation. As we dis- 
cuss below, the rate that the dipole moment relaxes to 
a constant is a measure of how the baryon density itself 
relaxes to its corresponding hydrodynamic modes. This 
rate will turn out to be much slower than the rate that 
the endpoints slow down, and this explains the discrep- 
ancy in the plot at times t ^ uj^/uq. 
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FIG. 11: A plot of the dipole moment of the baryon number 
and the string endpoint separation in the Minkowski spatial 
directions. The endpoint separation agrees very well with the 
dipole moment. There is however a small discrepancy near 
t ~ u\/uo, which in this plot is nTt ~ 17. This is due to the 
fact that the endpoint velocity slows down at a rate ~ e"'*'^^' 
whereas the velocity of the center of baryon density slows 
down at a rate ~ e"^'^"^'. 



Fig. 12 shows a plot of the angular distribution of 
baryon number 



dBs 
d6 



= 27r sin t 



^drps{t,x). 



(3.43) 



for a single jet as a function of both t and opening angle 
9. At distances close to the initial creation event dBg/dO 
is spread out. This just corresponds to the fact the lump 
of baryon density has a finite width and is very close 
to the origin. * In particular, it takes a time t ^ uq for 
light to reach the boundary after the initial event where 
the string is created. After a time i ~ a few mq, dis- 
tinct wave packets of baryon density corresponding to 
each quark will have formed and each packet will have a 
width ~ uq. From Fig. fTSl we see that after a time t ~ 



There is a potential numerical effect that might be influencing 
the spread of the ba ryon density at early times. When solving the 
equation of motion ( |3.9| for the radial component of the electric 
field E^{uj,q,u), we use a finite Fourier transform range with a 
cutoff A. Correspondingly in position space the numerically gen- 
erated baryon density cannot have a width that is parametrically 
smaller than 1/A. This effect can potentially make the early time 
behavior of dBs/dO artificially more spread out than it actually 
is. We have however checked t hat the early time behavior of 
dBs/d6 remains as shown in Fig.|12|as A is increased. 



a few M/i, dBs/dO has settled down to a relatively con- 
stant form. As discussed in Section IIA2| the time scale 
t ^ Ufi corresponds to the amount of time it takes for the 
string to evolve from the point-like initial conditions to 
the quasi-steady state configuration given in Eq. (2.251. 



During times in the window u/i <C i <C u\/uq, the baryon 
density slowly spreads out in the transverse directions. 
In Fig. [12] this manifests itself insofar that the angu- 
lar width of dBg/dO is approximately constant. Again, 
in this window of time the temperature dependence in 
the baryon density will come from that of the trajecto- 
ries of the string endpoints. That is, to leading order 
in uo/uh the angular width can be computed using the 
zero temperature expression for the angular distribution 
of baryon number in Eq. (3.361, but with the tempera- 
ture dependent endpoint trajectories. The precise result 
for the opening angle is rather complicated and depends 
on on details of the endpoint trajectory (i.e. the exact 
value of £,) so we do not give it here. 

At late times, the baryon density spreads out rapidly. 
In the gravitational theory this corresponds to times in 
which the string endpoints asymptotically approach the 
event horizon. As a string endpoint nears the horizon, to 
an observer on the boundary, the electric fields sourced by 
the string endpoint will have arbitrarily long wavelengths 
due to the strong gravitational redshift caused by the 
black hole. Correspondingly the baryon density induced 
on the boundary will be diffuse and spread out. From 
hydrodynamic considerations, one expects that the late 
time evolution of the baryon density will be that of the 
diffusion equation 



{do - DV^) p=Q 



(3.44) 



where D — 1/2ttT is the baryon density diffusion con- 
stant [2S1 [27] and is a localized source in space and 
time. Using techniques developed in Refs [5S1 |^ [50] . 
it is easy to compute the source from the gravitational 
theory. We have carried out this exercise, and for the 
symmetric strings discussed in Section |IIB[ found the 
following simple result 



e{t,x) = d-W6^{x)S{t) 



(3.45) 



where d is the dipole moment of the baryon density which 
as mentioned above, corresponds to the i — *■ oo limit of 
the string endpoint separation in the spatial direction. 
The above result for the source is easy to understand. 
At leading order in gradients, its form is fixed by the 
requirement that it be local in space and time and that 
the dipole moment of the baryon density in the hydrody- 
namic regime simply be d. 

We now turn to the rate in which dipole moment shown 
in Fig. [TT] relaxes to its final value of d. The hydrody- 
namic treatment of the baryon density discussed above 
always yields a constant dipole moment. Therefore, at 
late times, any time dependence in the dipole moment is 
a measure of how quickly the system relaxes to a hydro- 
dynamic description. Generically, one expects that the 
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FIG. 12: A plot of dBs/d9 for a typical string configuration where the endpoints travel far in the spatial directions. At times 
Uh ^ t <^ u\/uo, dBs/dO is highly localized about 6 = 0. During these times the baryon density wavepacket describes a 
quasi-particle. Around times t ~ u\/uo, which for this example is nTt ~ 17, the quasi-particle has lost most of its energy and 
the resulting jet thermalizes. Correspondingly, the baryon density of the jet spreads out and diffuses at late times. Note that 
the total baryon number contained in any slice of constant time is unity. 



rate that the system relaxes to hydrodynamics will gov- 
erned by the imaginary part of the quasinormal modes 
in the current-current correlator. For M = A SYM, these 
modes were first found in Ref |31j and are given by 



'2TTTn{±l + i), 



(3.46) 



The late time behavior of the center of charge can 
also be computed analytically from the bulk to bound- 



ary problem discussed in Section III We have carried out 



this analysis and found that the result simply reproduces 
the geometric optics argument given above. 



where n is a positive integer. We therefore expect the 
system to relax to hydrodynamics with a characteristic 
time scale l/27rT. 

As one can easily reason from the geodesic equation 
(2.211, the rate that the string endpoints slow down as 



they approach the event horizon is 



(3.47) 



so Xg{t) r^ e~ . This rate in fact applies under more 
general circumstances than for geodesies — any trajec- 
tory which is light-like will slow down at this rate near 
the horizon. Similarly, the rate that the radial coordinate 
increases is 



Us{t) - /. 



(3.48) 



Now suppose the source emits a signal at time ip- H will 
take some time t for this signal to reach the boundary. 
After a time At the source will have fallen a distance 



A5 



/ 



Ai 



At closer to the horizon. Therefore a 



signal emitted at time to -I- At will take a time r + 2 At to 
reach the boundary. It therefore follows that an observer 
on the boundary will see the string endpoint velocity in 
the spatial direction decrease like e~^'^^*. We therefore 
expect that the velocity of the center of baryon density 
decrease like e"^'^"^* instead of e"'*'^"^*. The time scale 
associated with the decay therefore agrees with the lowest 
quasinormal mode computed in Ref [31j 



IV. DISCUSSION 

A few remarks about our work are in order. We 
found that universally the endpoint motion of our strings 
rapdily settles down to a lightlike geodesic. Not just that, 
in fact the whole worldsheet approaches a lightlike con- 
figuration which we were able to construct analytically. 
There are two reasons why this was anticipated. For one, 
the effective string action goes to zero as the string falls. 
For the zero temperature string, the Nambu-Goto action 
in AdS is identical to Nambu-Goto in fiat space with a 
position dependent string tension that scales as ^. On 
top of this, the effective string tension also decreases as 
the string becomes more and more lightlike as it falls. 
For example, for both for a uniformly translating string 
X = vt aX zero temperature, and the dragging string of 
[HI |7] at finite temperature, the on-shell action scales as 
\\ — v^. One can again interpret this as a velocity de- 
pendent string tension that goes to zero as v goes to 1. 
Indeed, in our late time configuration each point on the 
string individually follows a lightlike geodesic. At finite 
temperature a similar point was also made in 0. 

One result we obtained that is perhaps surprising is 
that in general, the baryon density on the boundary is not 
localized directly above the string endpoints. In fact from 
Eq. (3.32 1 we see that at zero temperature and late times. 



the spatial velocity Vg of the string endpoints is different 
than the velocity Vs of the center of baryon density. The 
two velocities only agrees when Vs approaches the speed 
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of light. In this Hmit the baryon density on the boundary 
is localized directly above the string endpoints. 

To understand the origin of this behavior it is useful 
to understand the baryon density induced on the bound- 
ary by charges which exactly follow geodesies in the zero 
temperature AdS geometry. Consider two charges cre- 
ated at t = 0, which then move apart on the light-like 
geodesies 



Xs{t) = Vst 

u 



(t) = yr^Wt+ 



Uq. 



(4.1) 
(4.2) 



One can readily substitute the above trajectories into 
Eq. (3.24) and compute the induced boundary baryon 



density. A remarkable feature of the above geodesies is 
that the electric fields produced by the charges at times 
after the creation event do not induce any baryon density 
on the boundary — all of the baryon density is induced 
from the initial flash of light produced during the creation 
event. ^ Indeed it has been noted in several papers [ill^ 
that eternal geodesies, or geodesies which have infinite 
extent in the past and future, do not induce boundary 
densities at zero temperature. In the late time limit, the 
baryon density produced by the initial flash of light is 
simply 



pr(i,a;) = (-i) 



s+l 



(l-K^) 



Airr^il-x-Vs)^ 



S{t 



(4.3) 



This equation describes a focused spherical shell of 
baryon density expanding at the speed of light. While 
the shell of baryon density itself expands at the speed of 
light, the center of baryon density associated with each 
jet in fact moves at a speed less than the speed of light. 
This is because the baryon density is expanding in both 
the transverse and longitudinal directions. This is shown 
in Fig. 13 As one can easily compute from Eq. (4.3 1, 



the velocity of the center of baryon density is given by 
Eq. (3.321, and is in general not equal to Vs. 



At zero temperature, the fact that the baryon density 
induced by charges moving on geodesies is completely de- 
termined by the initial flash of light emitted during the 
creation event might lead one to suspect that the baryon 
density induced by string endpoints, whose trajectories 
only asymptotically approach geodesies, is very sensitive 
to the initial creation event and not just the asymptotic 
endpoint trajectories. This intuition is only partially cor- 
rect. While the value of the baryon density at a partic- 
ular point {t, x) might be sensitive to the details of the 
creation event and the histories of sources in the bulk, 
certain quantities which are averaged, such as moments 
of the baryon density, are not sensitive to the detailed 
history of sources in the bulk. One can explicitly see 
from Eq. (3.29) that the late time behavior of the center 



^ More precisely, all of the b aryon density comes from differenti- 
ating the e{t') in Eq. (|3.25i. 
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FIG. 13j A plot of the zero temperature baryon density given 
in Eq. (4.3 1 at three times ti, i2, is, with Vs = 0.8. At each 



time, the baryon density is localized on the solid black line. 
The coloring in the plot shows the angular distribution of 
baryon number, with red showing a high concentration and 
blue showing a low concentration. The yellow dots show the 
spatial position of the charge in the bulk at each time, with 
the top dot corresponding to time ts and the lowest dot cor- 
responding to time ti . The blue dots show the position of the 
center of baryon density at each time with the top dot corre- 
sponding to time ta and the lowest dot corresponding to time 
ii. As is evident from the plot, the position of the charge in 
the bulk does not coinci de with the center of baryon density. 
In fact, from Eq. (3.321 they move at different velocities so 



their separation grows linearly with time. 



of baryon density is completely determined by the late 
time behavior of the endpoint tr ajecto ries. Furthermore, 
one can explicitly see from Eq. (3.36) that the late time 



angular distribution of baryon number is also completely 
determined by the asymptotic behavior of the string end- 
points. Evidently, the information encoded in the radia- 
tion emitted during the times where the string endpoint 
trajectories are settling down to geodesies contains all 
of the information about the asymptotic behavior of the 
endpoint motion. 

For the case in which one string endpoint approaches 
a geodesic with Vs — 0, we see from Eq. (3.40 that the 



corresponding angular distribution of baryon number will 
be completely isotropic. Of course for a general value of 
Vs, one can always boost to a frame where this occurs. 
The angular coUimation of baryon density at late times 
is therefore a consequence of a Lorentz boost. 

We note that for the finite temperature string config- 
urations studied in this paper, the center of baryon den- 
sity did roughly coincide with the location of the string 
endpoint. We emphasize that the agreement is largely a 
consequence of the particular states we chose to study. 
We chose to study string configurations where the string 
endpoints traveled arbitrarily far in the spatial direc- 
tions before falling into the black hole. Such configura- 
tions naturally have the interpretation of quasi-particles 
in the boundary field theory. However, to generate such 
a configuration, the endpoint velocity must reside almost 
entirely in the spatial directions for most of the trajec- 
tory. Furthermore, the endpoints must move at the speed 
of light. Therefore Vs ^ I for most of the trajectory. 
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The fact that the baryon density resides over the string 
endpoint is therefore completely consistent with the zero 
temperature analysis given above. 

It would of course be interesting to study the stress- 
energy tensor of the jets produced by massless quarks. 
The stress-energy tensor of very massive quarks moving 
through a strongly coupled SYM plasma has been stud- 
ied in Refs [28 , 22J [321 El] while the stress-energy ten- 
sor of a meson moving through a strongly coupled SYM 
plasma has been studied in Ref |35j . However, while we 
have found that moments of the baryon density and the 
angular distribution of baryon number are in general in- 
sensitive to the exact initial conditions used to create the 
string, we suspect that the structure of the stress-energy 
tensor will be very sensitive to the initial conditions. The 
reason for this is that the strings corresponding to mass- 
less quarks in the field theory approach null strings after 
the initial creation event. As a consequence of this, the 
determinant of the worldsheet metric, which vanishes for 
a null string, is strongly dependent on the perturbations 
which exist on top of the null string. The 5d string stress 
tensor 



yields Ax 



v^r(f) 



1.31-^. Therefore a seemingly 



MN 



(Y) = - 



(4.4) 

will therefore also depend strongly on the perturbations, 
and consequently, via the gravitational bulk to boundary 
problem, the boundary stress tensor will as well. This 
should be contrasted with the electromagnetic bulk to 
boundary problem analyzed in this paper, where the elec- 
tromagnetic currents in the bulk J^ are not sensitive to 
the perturbations on top of the null strings. 

It would be advantageous to find averaged quantities 
related to the boundary stress-energy tensor, such as mo- 
ments or jet functions, which are not sensitive to initial 
conditions. If this is not possible, one can always in prin- 
ciple compute the boundary stress-energy tensor by in 
some sense averaging over a reasonable set of initial con- 
ditions or equivalently, by averaging over perturbations 
on top of a null string. Whether or not such an averaging 
procedure is manageable remains to be determined. 

At finite temperature one quantity which seems to 
be somewhat insensitive to the deviations from the null 
string is the total stopping distance Ax traveled by the 
quark. At least at the level of the null string approxima- 
tion, the stopping distance is determined by the initial 
radial coordinate uq and the geodesic parameter ^. For 
geodesies which only fall toward the horizon, the stop- 
ping distance can be obtained by integrating Eq. (2.231, 



Ax : 



du 

"0 vC 



1 



/(«) 



(4.5) 



This quantity is however, very sensitive to the value of 
^. For example, c onsid er the case of ^ = 1. Then in the 



Mo — > limit, Eq. (4.51 integrates to Ax = —. Now con- 



r(|) uo ■ «o ■ 
small change in ^ can yield a 30% change in Ax. The 
origi n of this behavior is simply that Ax as defined in 
Eq. (4.51, is not analytic at ^ = \fjijk))- Moreover, 
i — Vfi^o) sets the endpoint velocity in the spatial 
direction to be the local speed of light, and hence can 
maximize the total distance traveled. 

Based on the above arguments, one might suspect that 
the deviations from geodesies in the endpoint trajectories 
might have an important influence of the total distance 
traveled by the endpoint. From Fig. [7] one can see 
that the quantity, f{dxs/dt)~^, which is equal to ^ for a 
geodesic, is in fact not constant. It changes by 1 part in 
10^ over the course of the trajectory of the string end- 
point. This is consistent with the perturbative string so- 
lutions discussed in Section [II A 2[ In general, corrections 
to the endpoint trajectories will be 0{e). In fact, the ac- 
tual value that one chooses for ^ is arbitrary up 0{e) 
corrections. From Fig. |7] the possible values of <^ range 
from ^ - 1 = 6.5 X 10"^ to $ - 1 « -10"''. For the value 
of Uq — 0.06uh used in the numerical solution shown in 
Fig. [7] these two values for ^ lead to values for Ax which 
differ by a factor of two. Therefore, if one wants to ap- 
proximate the endpoint trajectory with a geodesic, one 
must carefully choose the value of £,. However, as shown 
in Fig. |7] ^ is essentially constant for the first half of the 
trajectory. It is during this part of the trajectory that 
the string endpoint is close to the boundary and the for- 
mula in Eq. (4.51 is most sensitive to the value of £,. If 
one uses the initial value of^ — 1 = 6.5x10"^, then one 
obtains a distance traveled by the geodesic which differs 
from the actual distance traveled by the string endpoint 
by 8%. 

Even if one chooses to use the initial value of ^ and 
approximate the trajectories of the string endpoints with 
geodesies, relating the parameters ^ and uq to the initial 
quark energy or momentum is in general nontrivial. This 
was discussed in Ref [5^, where the stopping distance of 
a gluon was considered. The analysis of [8, tried to re- 
late the parameters specifying the geodesic to the initial 
energy of the string. Without full string solutions only 
estimates were possible. Ref [8] estimated a stopping dis- 
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sider ^ = \/ f{uQ) and take uq — > 0. This configuration 



T^A ; """"" "- ~ ^ "■''^ ^string 

the energy of a string which asymptotically approaches 
the event horizon. In order to compare this estimate to 
our numerical strings, we take -^string to be half of the 
initial energy of our point-like strings. This is reasonable 
as the initial point-like configuration rapidly settles down 
to a quasi steady state configuration where the center of 
the string asymptotically approaches the event horizon. 
For the initial conditions given in Eq. (2.501, this yields 
-^'string ~ 1167\/Ar. Using the estimate of Ref [5], we 
therefore obtain a stopping distance of Ax « 13.3u;i. 
The actual distance traveled by the string endpoint is 
Ax « 19.3w/i. We therefore violate the estimate of Ref 
[B] by 30%. Consequently, we suspect that the analysis 
of Ref [SJ is incomplete. 
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Last, we discuss how our results generalize to string 
configurations which lie in multiple spatial directions. 
Such string configurations give rise to jets which are in 
general not back to back. Moreover, such configurations 
can give rise to additional jets caused by energetic gluons. 
As discussed in Ref [8^ , an energetic gluon is represented 
in gauge/string duality by a kink in a string. It is very 
easy to extend our analysis to general string configura- 
tions. Consider first the case of finite temperature. Sup- 
pose we consider a configuration which has in addition 
to the two jets caused by the string endpoints, N addi- 
tional gluon jets. Such configurations will have a very 
inhomogeneous distribution of momentum density in the 
initial string configuration. Portions of the string with a 
small momentum density will quickly fall to the horizon 
over a time scale t ^ u^- After this has happened, the 
trajectories of each portion of string extending up from 
the horizon are uncorrelated. The string endpoints and 
the gluon kinks will then move in arbitrary directions and 
each portion of the string profile will be very well approx- 
imated by steady state string profile. Therefore, the full 
string profile will look like 2 -I- A'^ trailing strings which 
are translating in different directions and are connected 
in a region of characteristic size ~ Uh- For times t^Uh, 
the baryon density induced by each string endpoint will 
therefore be virtually identical to that obtained by con- 
sidering strings which only exist in one spatial direction. 

To generalize our zero temperature results, we consider 
strings whose endpoint position in the spatial directions 
grows with time like \xg\ ~ t. Such strings will expand 
as they fall and will also approach a null string configu- 
ration. Moreover, the radial coordinate of each endpoint 
will increase with time like Us ~ t. We can therefore 
repeat the late-time expansion of Section |IIA 1| for this 
general case. Doing so, we find that the light-like bound- 
ary condition on each string endpoint forces its trajectory 
to asymptote to a light-like geodesic. Correspondingly, 
the angular distribution of baryon number for each jet 
will again be given by Eq. (3.36). However, there will in 



general be no correlation between the direction of each 
string endpoint and the resulting cones of baryon den- 
sity will thus point in arbitrary directions relative to each 
other. Moreover, the string asymptotes to a profile given 



by an arbitrary curve on the 4-sphere x 



= t\ So 



jets in large N^, strongly coupled Af — A super Yang- 
Mills theory at zero and finite temperature. Ultimately, 
we studied the baryon density of these jets. To do this, 
we obtained the endpoint motion for the dual strings. 
At zero temperature, a falling string's endpoint always 
travels a path that asymptotes to a lightlike geodesic. 
The finite temperature case is more complicated. Here 
we find that, for strings whose endpoints travel far in the 
Minkowski spatial directions, the endpoint motion ap- 
proximates a lightlike geodesic even at early times. We 
then numerically confirmed both of these results. 

Next , we computed properties of the baryon density for 
these jets. Our results are universal at zero temperature. 
In particular, two integrals of the baryon density — the 
center of charge of each jet and the angular distribution 
of baryon number — are completely determined by the 
endpoints' asymptotic trajectories. From this, we found 
that zero temperature jets travel forever and that their 
baryon density is focused into two cones whose angular 
dependence is related to the asymptotic motion of the 
string endpoints via the simple formula Eq. (3.401. 



as in the finite temperature case, a general falling string 
at zero temperature may have a number of kinks dual to 
gluon jets. 



At finite temperature, we see three time scales in the 
evolution of the baryon density. The first scale is the 
initial radial coordinate uq where the string was created. 
This scale sets the duration of time it takes light to prop- 
agate from the bulk of AdS to the boundary and sets the 
formation time for the production of quasiparticles. The 
second scale is the horizon radius u^. It takes a time 
t r^ Uh for the string in the bulk to relax to a quasi steady 
state configuration. The third time scale is uf^/uQ. For 
times in the window u/i <C t <C u\/uq, the string uni- 
formly translates along in the spatial directions and the 
endpoint slowly falls toward the black hole — the radial 
coordinate does not change much from its initial value of 
Mo. Correspondingly, during these times, the baryon den- 
sity is highly localized above the string endpoint and is 
naturally interpreted as the baryon density of a quasi par- 
ticle. At times t ^ u\/uo, the string endpoint is no longer 
close to the boundary and the resulting baryon density 
starts to spread out. This corresponds to the thermaliza- 
tion of the jet. For times t ^ u\/uq the string endpoint 
has fallen asymptotically close to the event horizon and 
the baryon density evolves according to hydrodynamics 
— it simply diffuses. 
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